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Abstract: We calculate the two-loop effective potential of the non-local NJL model de- 
rived from the Sakai-Sugimoto model in string theory. In contrast to conventional NJL with 
4-fermion contact interaction, the chiral symmetry was previously found to be dynamically 
broken for arbitrary weak coupling at the one-loop level. We calculate the one and two- 
loop contribution to the effective potential of the non-local NJL model and found that the 
two-loop contribution is negative. The two-loop potential for the chirally symmetric vac- 
uum is also negative but larger than the combined effective potential of the chirally broken 
vacuum at the two-loop level. The chiral symmetry breaking thus persists for arbitrary 
weak coupling. 
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1 Introduction 

Spontaneous Symmetry breaking (SSB) plays an important role in modern particle physics 
theory. Higgs mechanism in the standard model, for example, can be used to generate 
masses of elementary particles, leptons and quarks. The generation of quark masses by 
spontaneous symmetry breaking inevitably breaks the chiral symmetry of the QCD. Chiral 
symmetry breaking (xSB) could also be generated dynamically by the vacuum expectation 
value (vev) of chiral condensate ip{x)ip{y). Dynamical xSB can explain masses of mesons 
and hadrons which are responsible for most of the visible mass in the universe. It can also 
explain hadronic particle generation in strong interaction at low energies. The key idea 
of SSB is that any theory whose Lagrangian is invariant under some associated symmetry 
where vacuum state of such theory is not invariant and hence its vacuum carries non-trivial 
quantum number associated with the symmetry. In the spontaneous broken phase, there 
is an existence of Nambu-Goldstone (NG)-boson [1]. One can classify NG-boson into 2 
cases i.e. on one hand, NG-boson is massless elementary particle and on the other hand, 
such boson could be a composite particle. Dynamical symmetry breaking (DSB) usually 
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occurs as a result of the interaction between constituent particles in the theory and yields 
a composite NG boson. 

For chiral symmetry breaking of the QCD, the NG boson is usually identified with 
e.g. the three pions from the breaking of SU(2)l x SU(2)r to SU(2)v or the eight light 
mesons from the breaking of SU (3) l x SU(3)r to SU (3)y flavour diagonal. This symmetry 
breaking pattern [2] was successfully used to explore properties of the light hadrons and 
gives precise predictions of light hadronic spectra [3]. Early stage of xSB in the strong 
interaction was demonstrated by the linear-sigma model [4] and the current-algebra ap- 
proach [5, 6]. At the present, there is an incorporation between xSB and principle of 
effective field theory which gives a systematic framework to study QCD at low-energies, 
the so-called chiral perturbation theory [7]. The theory starts with an effective theory of 
hadrons with chiral symmetry in the action and use the SSB to generate a chiral symmetry 
breaking vacuum. The observed meson spectra shows good agreement with the prediction 
of the chiral perturbation theory [8, 9]. 

To address the chiral symmetry breaking/restoration phase transition, ones need to 
work with the action of quarks instead of hadrons. Nambu-Jona-Lasinio (NJL) model [10] 
is a model of quarks with four-fermion interaction which is employed to demonstrate the 
dynamical chiral symmetry breaking in the strong interaction independent of the confine- 
ment. Originally, NJL was formulated to explain mass of the nucleon as a consequence of 
the xSB. NJL model has been widely used as a description of low-energy effective model of 
hadrons in QCD with zero or finite-temperature [11-14] and electroweak symmetry break- 
ing by top-quark condensation or other fermions within or beyond the standard model 
[15-17]. NJL model is a very successful effective model to describe many hadronic proper- 
ties in low-energy QCD, for example, the mesons and baryons mass spectra, the pion decay 
constant, and the pion form factor (sec [11-13] for review). 

Despite the success of the NJL as a low-energy phenomenological model approach 
to low-energy QCD, the original NJL model does not address confinement. There are 
extensions of the NJL where inclusions of non-local interactions have been proposed in the 
literature (see [18] and references therein). One can simply reproduce the non-local NJL 
interaction from the QCD Lagrangian by integrating out the gluon field from the one-gluon 
exchanging diagram [15, 18]. 

In the non-local NJL approach, interaction depends on the momenta carried by the 
quarks leading to a momentum-dependent quark mass, generated by the spontaneous xSB. 
It has been shown that a non-local NJL model could lead to quark confinement with 
acceptable values of the parameters [19]. This phenomenon originates from the fact the 
quark propagator has no real poles and consequently quarks have no asymptotic states. 
There are several other advantages of the non-local NJL approach over the original (local) 
NJL model i.e. the nonlocality regularizes the model in a manner that anomalies [20] 
and gauge invariance [21] are preserved and the momentum-dependent regulator makes 
the theory finite to all orders in the 1/N C expansion. Finally the dynamical quark mass 
is momentum dependent in contrast to the original NJL model and consistent with lattice 
simulations of QCD [22]. As a result, one can see that the non-local NJL model may 
have more predictive power and be more realistic. There are two major applications of 
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the non-local NJL model in the strong interaction. Firstly, it is incorporated in the quark 
model to give mass spectra of excited mesons in good agreement with the experimental 
data [23]. Secondly, the thermodynamics of nuclear matter and QCD phase diagram could 
be explained quantitatively well by using non-local NJL model (with Polyakov-loop) [24]. 

A non-local NJL model can also be constructed from certain intersecting-branes con- 
figurations in string theory. The Sakai-Sugimoto model (SS) [25, 26] is a D8-D8-D4 
intersecting-branes model in type IIA string theory. The background spacetime is generated 
from a stack of N c D4-branes. An x A coordinate is compactified into a circle with radius 
R and the D4-branes wrap around the x 4 . On the boundary of the 10-dimensional space, 
a stack of Nf D8 and D8 are located at x 4 = —L/2 and L/2 respectively. The left (right)- 
handed quarks live on the D8 (D8)-D4 intersection in the form of open-string excitations. 
They are thus separated by distance L on the boundary and there is a U(Nf)L x U{Nf)R 
chiral symmetry. Geometrically, when the D8 and D8 merge at certain radial coordinate, 
the chiral symmetry breaking U(Nf)i, x U(Nf)n — > U{Nf)y occurs. 

We will not be considering the SS model in its full details in the present work but 
would rather focus on the low-energy effective 5-dimensional field theory limit of the model. 
In contrast to the strong coupling regime where the supergravity picture of intersecting 
branes provide simple geometrical interpretation of the theory, the weak coupling limit has 
its own unique picture of chiral symmetry breaking in terms of non-local NJL model in 5 
dimensions. 

In such intersecting branes setting there are two crucial parameters i.e. the 5-dimensional 
't Hooft coupling, A and the length scale of separation between D8-D8 flavor branes, L. 
One can consider the hierarchy of those parameter as A <C L which is the weak coupling 
regime. In such limit, we can treat left- and right-handed quarks as weakly interacting by 
single (five dimensional) gluon exchange process. The non-local NJL interaction is repro- 
duced by integrating out gluon fields in bulk spacetime from such D-branes configuration. 
In terms of effective potential in holographic non-local NJL, the nonzero solution of chiral 
quark condensate exists at arbitrary weak coupling [27, 28]. In contrast, if one considers 
the SS model in compactified case i.e. R is finite, and includes the KK tower of states. 
The xSB will happen only above a certain value of 't Hooft coupling [29]. In any cases, 
the analysis has been done on the effective potential of the non-local NJL at the one-loop 
level. It is interesting to investigate whether the two-loop contribution would change the 
profile of the effective potential in any significant way. 

We will start by reviewing the method of effective action in 5 dimensions when gauge 
fields propagate in 5 dimensions and fermions are localized in 4 dimensional subspace. By 
integrating out heavy gauge fields, we will obtain the effective fermionic action of the SS 
NJL model. Subsequently, by using auxiliary field approach, we integrate out the residual 
fermionic fields to obtain the effective scalar action of the SS NJL model. One-loop and 
two-loop contributions of the action are then calculated and discussed. Chiral symmetry 
breaking is demonstrated at both one and two-loop levels. 
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2 The effective Lagrangian 



We start with the effective action of the single-intersection model where left-handed quarks 
are located at a single intersection of N c D4 and Nf D8 branes [27], 



y = jd 5 x{-^F MN F MN + 5 (x 4 ) q\ <P (i + A„) 9L } , 



where M, N, 



0, 1, 2, 3, 4 and \i v , 



(2.1) 



0,1,2,3. Integrating by part and fix the 



gauge, the action can be rewritten in the following form, 

* = ^ / x { \ A m D AM + 5 jM A m) + / ^ x 1l i ^ 1L, (2.2) 

where we have defined = g\ q\ qi and set j( 4 ) = 0. 

The gauge fields live in 5 dimensions and it is natural to integrate them out to obtain 
4-dimensional effective action of the fermions. For consideration of the chiral symmetry 
breaking, we can bosonize the fermion bilinear and integrate out the fermions subsequently. 
In order to integrate out the heavy-gauge field Am- We recall the procedure from [30], start 
with 



jS J [dH] e * f I J [dH] 



,i / d 4 x3f (h(x),6) 



(2.3) 



where H{x) and l(x) are heavy and light fields respectively. In our case, Am is the heavy 
field and qL,R are the light fields. The actions with and without the light fields are given 
by 

j d 5 xJ?(A M ,q L ) = ^ J f 'd 5 x S^A M DA M + 8(x 4 )J M A M ^ 

+ J d A xq ] L ^id^q L , (2.4) 

/ d 5 x^(A M ,0) = \ [ d 5 x-A M aA M . (2.5) 
J 9% J 2 

By using functional path integral as demonstrated in appendix A, the effective action after 
integrating out the heavy gauge field can be read off from eq. (A. 3), 



y oS = iJ d^xql^d^qL 



16 vr 2 



d^xdSG(x-y^) q\(x)o»q L (y) q{(y) a„ q L (x) . (2.6) 



Next, we extend the Lagrangian (2.1) to the the left and right-handed quark fields located 
at different intersections, D4-D8 and D4-D8 respectively [27] , this is the low-energy field 
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theory limit of the SS model, 

y = f d5x { - ^ FmnFMN + 6 { x ' + ^) Ql^Wh + AJol 

= jj<? x {4? 9 mn □ 4 R) + * (* 4 + \ ) 4? + s (* 4 " |) 4? } 

+ y d 4 x q{ e»id^q L + J d A x q R a" i ^ q R (2.7) 

where we define JM\ = g\ q\ o M qi , = #5 gjj cr M q R and 4^ , are the gauge 

fields in 5-dimensional spacetime which are located on the D4-D8 and D4-D8 intersections 
respectively. 

Then the generating functional of the above action is given by 
f [dA$ d A$]A FP exp{iy{A%\A { K\q L ,q R )} 
= J[dA$>dAW]A FP 

xexp|^| d*xA^g MN UAf 

± J d 5 x d 5 y5 (x A + |) 6 (y 4 - |) jf L) {x) G MN (x - y , x A - y 4 ) jf R) (y) 
+ i J d A xq{a^id^q L + i J d*x q R a» ify q R J . (2.8) 

Using eq. (2.8), the effective action in the integrating out procedure is written by 

J[d A%> d A$ ] A FP expU^ (A$> , A$ ,q L ,q R )} 
f [d A$> d 4? ] A F p exp ji 5? (4? , 4? , 0, 0) } 

-ig£ I d A x d A yq ] L (x)^q L {x)G^{x-y,L) q R {y)a v q R {y) 1.(2.9) 



Finally, we obtain the effective non-local Lagrangian in the Feynman gauge as 
^cff = J d A x (q[ id fl q L + q R a* i <9 M q R ) 

-glj d A x d A y^g flu G(x-y,L)q l L (x)a fl q L (x)q R (y)a u q R (y), 
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d A x (ql a* id fl q L + q R ^ i q R ) 
+ A / ^ x d 4 yG(x-y,L)[q{{x)- q R {y)][q R {y)- q L {x)], (2.10) 

where we used the Fierz identity \cp L (x) <r M q^ixyj (q R (y) o"^ q R (y)J = — 2 (c^ L (x) ■ qniy)^ 

(qr{v) ' 9l( x ))- The dot in the right-hand side is the contraction in the colour indices, 
therefore each fermion bilinear in the final expression of the effective interaction Lagrangian 
is a colour singlet. There is a non-local interaction between two colour singlet operators in 
the theory. 



3 Effective potential at one-loop : Auxiliary field approach 

In this section, we will calculate the effective potential from the effective action eq. (2.10). 
Actually, there are many ways to perform this kind of calculation. We will use the standard 
method of effective field theory i.e. bosonize the fermion bilinear which would become the 
chiral condensate and integrate out the heavy-residual fields (in our case is the fermion 
fields) and then we can automatically obtain the effective potential with one-loop radiative 
correction from the effective Lagrangian. 

Following ref. [27], we start with the auxiliary field method. This method is used to 
study the symmetry breaking of the model by introducing the auxiliary field to the effective 
Lagrangian. In our case is the bosonized complex fields i.e. 

T{x,y) = ^-G{x-y,L)q\(x)- q R (y), 

T(y,x) = T^(x,y) = -^-G(x - y,L)q R (y) ■ q L (x), (3.1) 

where the coupling X/N c is related to the g\ coupling in the effective Lagrangian by the 

o 2 

relation A = -p? iV c . 

Substituting auxiliary fields from eq. (3.1) into the effective action eq. (2.10), we obtain 
y cS = J d A x (ql o» id^q L + q R a" i q R ) (3.2) 

+ J d 4 xd*y (~ X T G{x-yM + ^ X) ^ ' + T(g ' V) ^ ' qL{x) ) " 

In the chiral (Weyl) basis, one can rewrite the Lagrangian as 

J?W = J d 4 xq(x)(i$ + T(x)P L +T(x)P R )q(x)- J d'x^ ^J^ , (3.3) 

where we imposed the simplifying ansatz T(x, y) = T(\x — y\) consistent with the Poincare 
symmetry of the expectation value of the operator. This is justified since we are considering 
expectation value of T(x, y) in the vacuum to study the chiral symmetry breaking. 
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Figure 1. One- loop expansion of fcrmion fields. 



We are ready to integrate out the fermion fields in eq. (3.3), it reads 

jy eS _ f[ d qdq]e W {iy(q,q,T,T)} / J [dqdq] exp {i J?(q, q, 0, 0) } , 



exp < Tr In I 1 + 



T(x)P L + T(x)P R 



The identities /[dgdg] expji J d A xq{x)A q(x)} = det^4 = exp(Tr In .A) are used above. 
Then the effective potential with one-loop expansion can be determined from the effective 
action, 



Vr 



off 



off 



A * mm + iTr ln d + m ft + n*> ft , _ (3 . 5) 



A L) v * 

where Tr = Tr spinor Tr co i our Tr fl 

avor Tr S p ace ti me is the trace over all indices (i.e. spinor, color, 
flavor, spacetime). The physical meaning of this procedure is depicted by figure 1. 



The second term in the effective potential can be calculated by expansion 
t" m ) _ Tr S pj nor Tr co i or Trg avor Tr S p a , 



n=l 



"I) 



n-1 



n 



(TPl+TPr) 



iN c N f V I ^ilnfl + «« 



(2vr) 4 



k 2 



where we used the following relations; 



Tr S pi n or 3-spinor — 2 (in chiral basis) , Tr co i or l co ior — N c , 



Tr 



flavor -'-flavor 



Tl'spacetime = J d^X = 
74; 



j $ ifc- (z-j/) 

$ ~ J (2vr) 4 A; 2 



(3.6) 



(T#P L + T#P R y = TTk z . (3.7) 
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The momentum has been Euclideanized and henceforth we will drop the subscript. 
Finally, the effective potential at one- loop is given by (scaled by factor Nf) 



1— loop 



Nr. 



d 4 x T(x) T{x) 



.2" 



+ L 2 



A 



d 4 k 
(2 ttY 



In 



(3.8) 



The equation of motion (the gap equation) of the scalar T(x) from the effective action, 
eq. (3.8), is 



d x T(x) e 



-ikx 



(X 2 + P 

A 



T(k) 



k 2 + T(k) T(k) ' 



(3.9) 



Apart from the trivial solution T = for the chiral-symmetric vacuum, the general solution 
to the gap equation 8V e s/6T(k) = can be solved perturbatively either analytically or 
numerically (see appendix C). Non- vanishing T solution corresponds to chiral symmetry 
breaking vacuum which has lower energy and thus represents a true vacuum. We can 
obtain approximate solution by solving the gap equation in 2 regions of momentum, small 
and large k (i.e. T(k)T(k) S> k 2 and T(k)T(k) <C k 2 respectively). The two solutions then 
can be matched to determine the unknown constants. An approximate solution from such 
method is in the following form [27] 



T{k) = T{k) = < 



Tn = k* 



-<2 e" 



-o k 



< k < k* , 



S -S- , h < k < A. 



(3.10) 



Generically by using the gap equation, the one-loop potential can be rewritten to be 



1— loop 



N 



d 4 k E 



T(k)T(k) 
k 2 + T(k)T(k) 



(3-11) 



By substituting approximate propagator eq. (3.10) into eq. (3.11), we can demonstrate that 
there is chiral symmetry breaking vacuum induced by small momentum contribution to the 
one-loop potential. The details are discussed in section 5. Essentially, since the integrand 
in eq. (3.11) is a negative-definite function of variable k 2 /T(k)T(k), the one-loop potential 
is always negative for nonzero T regardless of the exact form of the solution of the gap 
equation. It is obvious that the solution with nonzero T gives the lower potential than the 
chiral symmetric solution T = 0. 

It is remarkable that the chiral symmetry breaking of the one-loop potential occurs at 
any weak coupling. The reasons are the boundness of the positive classical term (the first 
term in the right-hand side of eq. (3.11)) whilst the negative loop term (the second term in 
the right-hand side of eq. (3.11)) is not bounded for low momentum. The solution of the gap 
equation, eq. (3.10), is a constant for the low momentum, resulting in ln(l/fc 2 )-divergence 
of the loop term as k — > 0, regardless of A. 
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4 Effective potential at two-loop level 



Even though the one-loop potential demonstrates the possibility of chiral symmetry break- 
ing solution, higher loops contribution could likewise be significant. In this section, we 
will calculate the effective potential at two-loop level by following Jackiw's functional ef- 
fective action method [31]. The two-loop contribution can be calculated from the vacuum 
expectation value of the interaction Lagrangian Jfi(x) 



V; 



2— loop 



^OlTe^^-^lO), 



(4.1) 



where T is the time-ordering operator. In order to obtain the two-loop contribution, we 
simply use the conventional Feynman rules to calculate all possible two-loop diagrams exist 
in the effective theory of fermion and auxiliary scalar. The propagator Q of the field (f) to 
be used in the evaluation of the 2-loop diagrams is defined by the inverse of the functional 
operator, namely 



iQ 1 (x,y) 



8(/)(x) 6c/)(y) 

In the previous section, we recall the effective Lagrangian 



(4.2) 



off 



<?X 



i q(x) $ q(x) + q(x) T(x) P L q(x) + q(x) T(x) Pr qix) 



N c T{x) f(x) 
A G (x, L) 



5?! = / d 4 xJ^/(z) 



d A x 



q(x) T(x) Pl q(x) + q{x) T(x) Pr q(x) 



(4.3) 



The functional operator S 1 of the quark fields and T> 1 of the complex scalar fields from 
the effective Lagrangian therefore can be written as 



iS-\k) = J d 4 xe- ik < x ~ y) 
= $ + T{k) P L + f 
iV~ l {k) = [ d 4 xe~ ik < x -^ 



6q(x) 5q{y) ' 
+ T(k) P L + T(k) P R , 



5T(x)ST(yY 



N r . 



XG(k, L) 



(4.4) 



Noting that, T> 1 has no kinetic term for the scalars T,T . The vertices of the interaction 
in the effective action are given by [31, 32] 



Tr - [ J* rr -ik-(x-y-z) (-ifS 3 y cS 

J 6q(x)5q(y)T(zy 
= iP L , 

r R = f d 4 xc- ik - (x - v ~ z) (~^ 3 / 3 ^ ff 

J Sq(x)Sq(y)T(z)' 
= iPB ■ 



(4.5) 
(4.6) 
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Putting everything together, the two-loop contribution is (see figure 2 for the corresponding 
diagram) 

d A k d 4 p 



V 2 -ioo P - *Tr , 4 4 



_A_ 



r L s(p)T R s(k)v( P -k), 

d A k d 4 p 2p ■ k G (p - k, L) 



NfNc 1 (2^) 4 (2vr) 4 [p 2 -T(p)f(p)] [k 2 -T{k)f(k)\ 



where we have used 



ft-TP L -TP R 



(f + TP L + TPr) (f-TP L - TPr) 



TP L -TP R 
k 2 -TT 



and 



Tr | Pj, (p 1 — T Pl — T Pr) Pr (ft — T Pl — T Pr) j = 2 



p • A;. 



(4.7) 



(4.8) 



(4.9) 



Using Wick rotation i.e. k ik E , d A k — > i d 4 kE , we obtain the two-loop contribution in 
Euclidean space 



d A k E 



1 



V2-ioo P = -4^XN f , (27r)4 

d A p E 2 PE -k E e- L \ PE ~ kE 



k 2 E + T(k E )f(k E ) 



(2 7T) 



p 2 E + T{p E )T{p E ) 



Pe ~ k E 



(4.10) 



The angle integration can be evaluated as shown in appendix D to be 



V; 



XNt 



2-loop - 4L57r 3 



LA 



dkE 



k E 



k 2 E + L 2 T(k E )f{k E ) 



LA 



dpi 



- i 
Pe 



p 2 E + L 2 T(p E )f(p E ) 



E 

n=0 



(-l) n I ir(n-l)A 



a — .-j 
2 



(n + 3) (n + 5)B 



2 (A 2 - B 2 ) 2 Fi 



i 3 - n 5 - n B 2 
4 ' 4 ' A 2 



/ ^9 / \ ,o\ _ / 3 — n 5 — n -B 2 
+ (i? 2 (n + 2)-2A 2 ) air^—,—^;-^ 



(4.11) 



where kE = L kE , Pe = L pe , A = p E + k E , B = 2pE kE • Henceforth, for convenience 
we will simply write the Euclidean momentum without a subscript. 

5 Results and Discussions 

Adding all of the one-loop and two-loop contributions, the total effective potential becomes 

—loop 

+ v 2 

— loopi 



Nr. 



d A k E 
(2tt)4 



T(k)f(k) T(k)f(k) 



k 2 + T(k)T{k) 



k 2 



2— loop • 



(5.1) 
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Figure 2. Two-loop vacuum diagram for V2-I00P, 
scalar. 



solid line is the fcrmion and dash line is the 



The resulting effective potential of the scalar shows the possibility of chiral symmetry 
breaking at nontrivial T{k) 7^ since the sign of the one-loop contribution is opposite to 
the classical action of the scalar. The 2-loop effect as given in the form of eq. (4.11) could 
be either positive or negative depending on the relative sizes of each n-term. A closer 
investigation reveals that the n = term is the largest and it is negative. The odd-n 
terms are positive with smaller values than the preceding even-n terms. Consequently, 
the entire two-loop potential is negative. Since the chiral symmetric solution T = T = 
gives larger negative two-loop contribution than the chirally broken case (with smaller 
denominator of the integrand in eq. (4.11)). It is thus possible that the difference of 2-loop 
contributions would compensate the one-loop effect and alter the true vacuum of the theory 
in a significant way. We will demonstrate that the two-loop contribution is small comparing 
to the leading one-loop and the chiral symmetry breaking persists. In evaluation of the 
momentum integrals, we will apply a UV-cutoff A required in non-renormalizable effective 
field theory. The cutoff will be taken to be larger than To and smaller than 1/L. 

5.1 1-loop 

Since both one- loop and two- loop contributions scale with the number of flavour Nf, we 
will simply suppress the Nf factor henceforth. First, we will consider 1-loop contribution 
of the scalar to the effective potential and demonstrate that the potential has nontrivial 
minima when using the ansatz solution of the gap equation as given in eq. (3.10). The 
1-loop integrations, eq. (3.11), can be performed in two separate momentum regions and 
rewritten as the following 



To 

+ 



T{k)f{k) 
k 2 + T(k)T(k) 



In 1 + 



T 



(k)f(k) \ 
k 2 







T 



d A k 

(2tt)< 



T(k)f(k) 
k 2 + T(k)T(k) 



In 1 + 



T{k)f{k) 
k 2 
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rpi 
J 



In 2 



16vr 2 v 

where we define the function 



1 

n=l 



\n-l 



(1 



F 4n _ 3 (LT ,LA;n) 



(5.2) 



F 4n _ 3 (LT ,LA;n) = E 4n _ 3 (2nLT ) - E 4n _ 3 (2nLA) 



4rt-4 



E m {z) 



oo — 2t 



n & 

The function i ? 4 n _3(LTo, LA; n) decreases very rapidly with n, therefore the sum in the 
one-loop potential, eq. (5.2), can be approximated by truncating at finite n with a high 
precision. 

The one-loop contribution will be explored by fixing one and two of the 3 param- 
eters, A, L, A and numerically plot the effective potential with respect to the remaining 
parameters. The physically- valid region of the parameter space for our SS NJL model is 
To = a/A/L 3 < A < 1/L. As shown in figure 3, the one-loop potential is negative at any 
nonzero values of A, A, L corresponding to nonzero values of To = \J A/L 3 . Since when 
T = 0, the potential is zero and less preferred than negative potential occuring at any 
coupling A, chiral symmetry breaking thus naturally occurs for any weak coupling (i.e. 
To < 1/T). In figure 3(a), the potential approaches negative constant for A > 0.7 — 0.8 
for A = 0.1, L = 1. If we instead fix the UV-cutoff scale A = 0.5, the potential will be 
a decreasing function with A as demonstrated in figure 3(b). Figure 3(c) also shows the 
one-loop potential at fixed A = 0.3, A = 0.5 as a function of L. It is important that we 
restrict ourselves to the physical region To < A < 1/L in our consideration of the effective 
potential in the nonlocal NJL model. 

We also plot the potential landscape in the physical region at fixed L = 1, as is shown 
in figure 4. 

5.2 2-loop 

In this section, we investigate the 2-loop contribution to the effective potential. Using the 
one-loop approximate scalar ansatz, eq. (3.10), the two loop integration given by eq. (4.11) 
can be separated into 3 terms, 

A 



2— loop 



4vr 3 L 5 



T L _ rT L r-AL _ ML pT L _ f-AL 

dk dp+ dk dp +2 dk dp 

Jo JtqL JtqL Jo Jt q l 



p 



\ k 2 + L 2 T(k/L) T(k/L) ] [p 2 + L 2 T{p/L) T(p/L) 
~ (_l)n / 7r(n-l)A^ 



n=0 



(n + 3) (n + 5) B 



2 {A 2 - B 2 ) 2 Fi 



3 — n 5 — n _ B A 



4 " ;1; A 2 



I 3-n 5-n B 2 
4 ' 4 ' 'A 2 



+ (B 2 (n + 2)-2A 2 ) 2 Fi 



(5.3) 
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Effective potential versus UV -cutoff Effective potential versus A 




1.1 1.2 1.3 1.4 1.5 
L 



(c) 

Figure 3. One-loop effective potential per colour as a function of A, A, L. 




Figure 4. One-loop effective potential per colour as a function of A, A for L = 1. 

The overall 2- loop contribution scales with A/L 5 . The integration in the low momentum 
region has additional (TqL) g x (ToL) 71 " 1 factor for each n-term in the sum. The integration 
in the high momentum region, on the other hand, has additional (AL) 6 x (AL) n_1 de- 
pendence for each n-term. The cross term integration has additional overall scaling factor 
(ToL) 3 (AL) 3 for all n. We perform numerical integration on each n-term and add them 
up. Since the integrand for each n is a smooth and well-behave function with no singulari- 
ties and abrupt changes, numerical integration yield very precise results. The value of the 
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integration for each n-term decreases rapidly with n and the error is less than 10 -6 if we 
truncate the sum at n = 10. 

Figure 5 shows the effect of 2-loop contribution to the effective potential. Similar to 
the one-loop case, the chiral-symmetry broken vacuum solution has lower energy than the 
chiral symmetric one (T = T = 0). However, there is one crucial difference between one and 
two-loop potential. The magnitude of 2-loop contribution could increase with the cutoff 
in contrast to the 1-loop which saturates to negative constant. This is originated from the 
AL dependence of the potential in eq. (5.3) getting larger with increasing A resulting in 
the decreasing function of the effective potential with the cutoff. 



Effective potential versus UV-cutoff 



■S -0.00005 - 



-0.00006 



N c = 3, A = 0.1, L= 1 



0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Or 

-5.X10- 7 " 

-1.x 10"' 
-1.5x10-' 

-2.x 10-' 
-2.5X10" 6 

-3.X10"'- 



Effective potential versus UV-cutoff 



-loop ^2-loop 



iV c = 3,A. = 0.01, L= 1 



0.2 



(a) 



(b) 



Potentials versus coupling A for L=l, A=0.9 




(c) 

Figure 5. Effective potential up to 2-loop as a function of A for N c = 3, A = 0.1 (a), 0.01 (b); 
L = 1 . The 2-loop contributions for both chirally broken and symmetric solutions are shown in (c) 
in comparison to the chirally broken 1-loop. 



When the 't Hooft coupling is very small, A < 0.08, interesting phenomenon occurs. As 
we can see from figure 6 (a), the ratio (magnitude) of the 2-loop to 1-loop increases sharply 
as A — > 0. From eq. (5.2) and (5.3), since Vi-ioop ~ A 2 ^ while V2-I00P ~ A, the ratio of the 
2-loop to 1-loop potential will scale as (XN c )~ l and the 2-loop contribution will be dominant 
for sufficiently small A. However, since both one and two loop contributions are negative 
and together they are larger in magnitude than the two-loop potential of the chirally sym- 
metric solution (figure 5 (c)), the chirally broken vacuum always has lower potential as is 
demonstrated in Figure 5 (c) and 6 (b). Chiral symmetry breaking therefore persists for 
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0.00 0.01 0.02 0.03 0.04 ' 0.00 0.02 0.04 0.06 0.08 0.10 



A A 
(a) (b) 

Figure 6. The ratio of V^-ioop/^i-ioop decreases rapidly with increasing A (a). In (b), the effec- 
tive potential for small A remains negative (and smaller than the chirally symmetric case), chiral 
symmetry breaking thus persists. The cutoff A is set to 0.9 and L = 1 for these plots. 

arbitrary weak coupling. One might anticipate the chirally symmetric potential to become 
more negative than the chirally broken one as A — > since V2_i 00 p(sym) < V2_i 00 p(xSB) 
and V2-ioop/^i-ioop ~ (-^cA)" 1 . However, a closer investigation reveals that the differ- 
ence V^-ioopCsym) - V 2 _ioo P (xSB) ~ AV 2 -ioo P (sym) for very small A and Vi-\ oop (xSB) is 
actually larger in magnitude (i.e. more negative) than AV2-ioop(sym). Consequently, the 
chirally broken solution still has lower potential than the chirally symmetric one even for 
extremely small coupling. 

6 Conclusions 

SS intersecting-branes model provides a geometrized model of chiral symmetry breaking 
and confinement in both weak and strong coupling regimes. The effective field theory at 
low energy (E < 1/L) from the SS model is a type of NJL model with non-local 4-fermion 
interaction. In constrast to conventional NJL with 4-fermion contact interaction which 
requires sufficiently large coupling to break chiral symmetry, the holographic non-local 
NJL model prefers chirally broken phase for arbitrarily weak coupling at the 1-loop level. 
In this work, we found that the 2-loop effect does NOT change this feature of the model 
and the chiral symmetry breaking persists for arbitrary weak coupling. The 2-loop effect 
can be understood as the antiscreening of the non-local 4-fermion interaction induced from 
the coupling of fermion with the cloud of scalar condensate. The 2-loop contribution in our 
model is suppressed below the 1-loop contribution for A > 0.08. One of the suppression 
factor is the number of colour degrees of freedom N c = 3 when the coupling A is fixed. 
Bosonization of the fermion bilinear into a colour-singlet scalar naturally matches the loop 
expansion of the effective potential with the 1/N C expansion. The large- N c expansion 
makes it manifest that the 1-loop contribution scales as (Nc) 1 and the 2-loop scales as 
(N c )°. The higher-order loops are therefore suppressed by negative power of iV c and so on. 

The solution to the gap equation in our model gives fermion condensate proportional 
to \/X/L 3 resulting in A 2 -dependence of the one-loop potential while the two-loop scales as 
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A. Therefore the two- loop contribution becomes dominant to the one-loop for very small 
coupling. However, since the difference between the two-loop contribution of the chirally 
symmetric and broken solutions is smaller than the size of the one-loop potential of the 
chirally broken solution, chiral symmetry remains broken for arbitrary small coupling. 
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A Integrating out heavy gauge field for single intersection model 

According to the standard technique of path integral (see [30], or textbooks in QFT), the 
generating function of the action in eq. (2.4) is written by 

[dA M ] A FP exp{iS'(A M ,q L )}, 

J[dA M ]A FP exp {z j d 5 x^ i A M □ A M + 5 (x 4 ) J M A M + i J d 4 x q\ a^id^q^, 



j [dA M ]A FP 

x exp { 4" / d 5 x - A M □ A 
1 91 J 2 



^ / d 5 x d 5 y5(x 4 )5(y A )\j M (x)G M N(x-y,x A -y*)J N (y) 

95 J 1 



+ ijd 4 xq[ a»id^q L y (A.l) 



where A^p is the Faddeev-Popov's determinant. In the Feynman-gauge, the propagator, 
Gmn(x,x 4 ) of the gauge field Am can be written as 



G M n (x , x 4 ) = — r g M N G (x , x 4 ) , 

1 gMN 



2 (tryi\2 _ ^.2^/2- 



l7T^ (( x 4)2 _ x 



(A.2) 
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Using eq. (2.5,A.l) in eq. (2.3), we obtain 

J*<* = J [dA M ]A FP e^ d4 ^( A ^^) / J [dA M ]A FP e^ dix ^ A ^°) 
[dA M ]A FP 



\ 4 / d 5 x^A M OA M 



9l 



\ I d 5 xd 5 yS(x 4 )S(y 4 )±J M (x)G M N(x-y,x 4 -y 4 )J N (y) 
91 



+ i J d A xq\ a^id^q^ 



[dA M ] A FP exp , 

95 



d 5 x - A M □ A 



}■ 



ex P { ~ ~2 J d ^ x Yg^2 95 Ql( x ) Ql(x) G(x -y,0)gl q[(y) q L (y) 

+ i J d 4 xq j L u^id^qi^, 



exp 



{-ijd^x d% J^G(x-y,0) [q[(x) a* q L (y)] [q[(y) <r„ q L (, 



d*xql*' l id ll q L }, (A.3) 
where we used J (4) = and Fierz identity (q* L>1 q Lt2 ) (q^ 3 ^ q LA ) = (q\ x q LA ) (q[ 3 q L:2 ) 

B Fourier transform in Euclidean 5-dimensions 

The coordinates of a d dimensional Euclidean space are given by 

x\ = x cos 6\ , 

X2 = x sin #i cos 02 , 

X3 = x sin 9\ sin 02 cos 0s , 

£4 = x sin 6*i sin 02 sin #3 cos 6*4 , 

X5 = x sin 6*i sin 02 sin #3 sin #4 cos #5, 



x n = x sin 0i sin 02 sin #3 • • • sin0 n _i cos0 n 
x n+ \ = x sin$i sin 02 sin #3 • • • sin0„_i sin0 ra , 



(B.l) 



where only n G [0, 2ir) (so that x n +i S [— x, +x, ]) and other angles range from to ir. For 
Euclidean momentum in 5 dimensions, the components are 

k\ = k cos 0\ , 

k2 = k sin 0\ cos 02 , 
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k 3 = k sin 0i sin 2 cos 9 3 , 

k A = k sin 61 sin 02 sin #3 cos 84 , 

k§ = k sin#i sin 9 2 sin #3 sin #4. (B.2) 

The volume element in 5-dimension is thus 

d 5 k = k 4 sin 3 0i sin 2 # 2 sin <9 3 dfc c26>i d6 2 d0 3 d9 A . (B.3) 

With this measure, the Fourier transform in 5 dimensional Euclidean space can be 
performed as the following 

F{x) = I iw p{k)elk ' x= I iw p{k)elkxcosei > 



1 



(2vr 



1 5 



dk dOi d0 2 d6 3 dO A k 4 sin 3 61 sin 2 9 2 sin 9 3 F(k) j k x cos 01 , 



(27T) 

1 

47T 3 



/"OO /"TT /"TT /"IT r2lT 

— F / F(k)k 4 dk / e ife:ECOs91 sin 3 ^^ / sin 2 2 <W 2 / sin# 3 d# 3 / d0 4 , 



Akx cos 6* 



Given an explicit functional form F(k), the transform can be completed. 

However, in the situation where the quarks are localized at particular x 4 and the gauge 
fields in 5 dimensions are integrated out to obtain the effective 4-dimensional action, we 
will need to perform the Fourier transform of the propagator given in eq. (A. 2) under 
the condition that the gauge fields are propagating at a fixed distance in x 4 direction. 
The Fourier integration will split into a delta function in x 4 coordinate and the Fourier 
transform in the Euclidean 4 dimensions. 

For example, in our model, the Fourier transform becomes 

/poo r 
d 5 xG(x,L)e ik - x = / dx A e ikiXi I d 4 x G(x , L) e~ ik ' x , 
J — OO J 

= 5{k A ) j dx d9\ d6 2 d6 3 x 3 sin 2 #1 sin 6 2 G(x , L) e 

= 6 (k A ) w x 3 dx 7rJl } kx) (2) (2tt), 

Jo (L2 + 5 2)2 kx 

47T 2 [°° {kxfj^kx) 

k Jo {{kLf + {kx) 2 Y 

— Lk 

= 5{k A )AK 2e —, (B.5) 

where J n (x) , is Bessel function. If we neglect the momentum in bulk spacetime say, k A = 0, 
we obtain 

— L k 

G(k,0) = 47T 2 ^— . (B.6) 
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C Gap equation at two-loop 



In this section, we will derive the gap equation at the two-loop level. Start with the 
two-loop effective potential 



V eS = N, 
- A 



d i xT{x) T(x] 



(x 2 + L 2 ) 



d 4 k 
(2 ttY 



In 1 + 



T(k)T(k) \ 
k 2 ) 



d 4 k 



(2vr) 4 J (2tt) 4 



d A p G(p - k, L) 



2p • k 



k 2 + T(k) T{k) p 2 + T{p) T(p) 



(C.l) 



The functional derivative of V e s with respect to T(q) gives the gap equation 



SV e 



cff 
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ST(q) ST(q) 



N c / d A xT(x) 



d 4 k 



d 4 k 



In 1 + 



-ik-x 



f(k) 



d 4 k 



(2vr) 4 J (2tt) 



T(k)T(k) \ 
k 2 J 

d 4 p G(p - k, L) 



2p-k 



k 2 + T(k) T{k) p 2 + T(p) Tip) 



f J d 4 xTix) J ^e-*-*6W(k- q )(x 2 + L 2 )l 



T(k) 



d 4 k 

(2tt) 4 k 2 + T(k) T(k) 
d 4 k f d 4 p 



S^(k-q) 



2p-kGip-k, L) 



(2vr) 4 J (2vr) 4 

k 2 Tjp) S^jp-q) 

[k 2 + T(k) f(k)] 2 [p 2 + Tip) f(p)] 2 [k 2 + T(k) f(kj\ 2 [p 2 + T{p) f(p)] ' 



T(k)T(k)T(p) 5^ip 



Tip)Tip)Tik)5 { - 4 \k 



p- 



Tik)5^ 4 \k-q) 



[k 2 + T(k) f(k)] 2 [p 2 + Tip) Tip)] 2 [k 2 + T(k) f(k)] 2 [p 2 + Tip) Tip)] 2 J ' 



Nr. 



A(2vr) 4 
A 



+ 2- 



d 4 x e~ iq - x 
d 4 k 



N r 



Tig) 



Tjx) 

G(x,L) (2vr) 4 <?+T(q)T(q) 
2q-kGiq-k, L)k 2 Tiq) 



+ 2- 



(2tt) 4 7 (2^) 4 [h* +Tik)T(k)] 2 [q 2 + Tiq)fiq)]' 
d 4 k 2q ■ k Giq - k, L) T(k) T(k) T(q) 



(2 7r) 4 7 (2vr) 4 [ fc 2 + T ( fc ) f (fc)] 2 [ g 2 + T ( ? ) T ( ? )] 



2 ' 



(C.2) 



where we have used the symmetric property of the propagator L) i.e. — p, L) 
G(\k-p\,L) =G(]p-k\,L). One then obtains the E.O.M. iSV e g/df(q) = 0) as 



Gix,L) 



N r 



q 2 + T(g) T(g) 



+ 2A 



d 4 k 



2q ■ kG(q - k, L) 



k 2 Tiq) 



[k 2 + T(k)Tik)] 2 [ q 2 +Tiq)T(q)] 2 
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+ 2A [ d ' k 2a-kG(a k L) nk)T{k)T{g) 
C.l The /c 2 > T(k)f(k) approximation 

We will approximate the gap equation in two regimes. First when k 2 3> T{k) T(k), eq. (C.3) 
becomes 



G(x,L) q 2 

+ "/(0 29 '* G( '-*' i) W i 

^/^.^-^TOm^. (C4) 

Multiply by and neglect the last term in the left-hand side, we obtain 



A 



/ ^ ^ ^Jl) " ^ q2T{q) +2X j 2 <l-kG(q-k,L)^- = 0(.C5) 



The last term represents the non-local screening effect of the scalar which is ^-suppressed 
comparing to the other terms. If we neglect the screening effect and Fourier transform the 
rest, the one- loop gap equation is recovered, 

v2 (<|B)) +ATW = ' <c - 6) 

where V 2 is the Euclidean Laplacian in 4 dimensions. 
C.2 The T(k) f(k) > k 2 approximation 

Next, we consider to the low momentum regime i.e. T{k) T(k) 3> k 2 , the E.O.M. in this 
limit is given by 

Nc f fxe-*-* ?} X \~Nc T[q) 



A7 G(x,L) "T{q)T{q) 

d'k n _ t r , k 2 T(q) 



+ 2A j ■^2q.kG(q-k, L) 



[T(k)T(k)} 2 [T(q)f(q)] 2 



+ 2x[**2 q .kG( q -k,L) M 2 = 0. (C.7) 

J ^ [T(k)T(k)) 2 [T(q)T(q)] 2 

Neglecting the third term in the left-hand side, the gap equation becomes 

X J G(x, L) ~ iVc f(q) J (2 vr) 4 T(k) f(k) T(q) f(q) f(q) ' ^ 

The last term on the right-hand side represents the non-local screening effect of the scalar 
which is Ar c -suppressed. Again, the one-loop gap equation is recovered when the screening 
effect is neglected. 

The gap equation at one-loop level is solved in ref. [27] as given in eq. (3.10). We use 
this approximate solution in the evaluation of the effective potential. 



- 20 - 



D Evaluation of the 2-loop angle integration 



We will integrate out the internal angle of Euclidean 4-dimension, we start with the two- 
loop effective potential, 



^2-loop = -47T 2 AiV) 



Ski 



(2vr) 4 
4vr 2 AiV- rA 
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dk E t 



k 2 E + T(k E )f(k E ) 
(2vr 2 )4 



J (27T 



Pe 2pe • k E 
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where k E = L k E , p E = L p E , A = p 2 E + k E , B = 2p E k E 

Expanding function ^ = \ Yln=o 7T\ § ives 
\N- rLA 
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OG 



(-l) n ir(n-l)A 

x L 



re — 3 
2 



n ! In + 3) (n + 5) B 

, .o „o, „ / 3 — n 5 — n i? 2 

2(^ 2 -5 2 ) — ;i;^ 

+ (S'(n + 2)-2A>) ^(^,^2;^ 
where 2-^1 («> &j c > -2) is the hypergeometric function. 
References 



(D.2) 



[5] 
[6] 

[7 



[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 



J. Goldstone, Field Theories with Superconductor Solutions, Nuovo Cim. 19, 154 (1961). 

C. Vafa and E. Witten, Restrictions on Symmetry Breaking in Vector-Like Gauge Theories, 
Nucl. Phys. B 234, 173 (1984). 

J. Beringer et al. [Particle Data Group Collaboration], Review of Particle Physics (RPP), 
Phys. Rev. D 86, 010001 (2012). 

M. Gell-Mann and MLevy, The axial vector current in beta decay, Nuovo Cim. 16, 705 
(1960). 

M. L. Goldbcrger and S. B. Treiman, Decay of the pi meson, Phys. Rev. 110, 1178 (1958). 

M. Gell-Mann, R. J. Oakes and B. Renner, Behavior of current divergences under SU(3) x 
SU(3), Phys. Rev. 175, 2195 (1968). 

J. Gasser and H. Leutwyler, Chiral Perturbation Theory to One Loop, Annals Phys. 158, 142 
(1984). 

S. Scherer, Introduction to chiral perturbation theory, Adv. Nucl. Phys. 27, 277 (2003) 
[hcp-ph/0210398]. 

H. Leutwyler, On the foundations of chiral perturbation theory, Annals Phys. 235, 165 (1994) 
[hep-ph/9311274]. 

Y. Nambu and G. Jona-Lasinio, Dynamical Model of Elementary Particles Based on an 
Analogy with Superconductivity. 1., Phys. Rev. 122, 345 (1961). 

U. Vogl and W. Weise, The Nambu and Jona Lasinio model: Rs implications for hadrons 
and nuclei, Prog. Part. Nucl. Phys. 27, 195 (1991). 

S. P. Klevansky, The Nambu- Jona-Lasinio model of quantum chromodynamics, Rev. Mod. 
Phys. 64, 649 (1992). 

T. Hatsuda and T. Kunihiro, QCD phenomenology based on a chiral effective Lagrangian, 
Phys. Rept. 247, 221 (1994) [hcp-ph/9401310]. 

M. Buballa, NJL model analysis of quark matter at large density, Phys. Rept. 407, 205 
(2005) [hcp-ph/0402234]. 

V. A. Miransky, Dynamical symmetry breaking in quantum field theories, World Scientific, 
Singapore (1993). 

W. A. Bardccn, C. T. Hill and M. Lindner, Minimal Dynamical Symmetry Breaking of the 
Standard Model, Phys. Rev. D 41, 1647 (1990). 



- 22 - 



[17] C. T. Hill and E. H. Simmons, Strong dynamics and electroweak symmetry breaking, Phys. 
Kept. 381, 235 (2003) [Erratum-ibid. 390, 553 (2004)] [hep-ph/0203079]. 

[18] G. Ripka, Quarks bound by chiral fields, Oxford University Press, (1997). 

[19] R. D. Bowler and M. C. Birsc, A Nonlocal, covariant generalization of the NJL model, Nucl. 
Phys. A 582, 655 (1995) [hep-ph/9407336]. 

[20] E. Ruiz Arriola and L. L. Salcedo, Chiral and scale anomalies of nonlocal Dirac operators, 
Phys. Lett. B 450 (1999) 225 [hep-th/9811073]. 

[21] B. Golli, W. Broniowski and G. Ripka, Solitons in a chiral quark model with nonlocal 
interactions, Phys. Lett. B 437, 24 (1998) [hep-ph/9807261]. 

[22] M. B. Parappilly, P. O. Bowman, U. M. Heller, D. B. Leinweber, A. G. Williams and 

J. B. Zhang, Scaling behavior of quark propagator in full QCD, Phys. Rev. D 73, 054504 
(2006) [hep-lat/0511007]. 

[23] R. S. Plant and M. C. Birse, Meson properties in an extended nonlocal NJL model, Nucl. 
Phys. A 628 (1998) 607 [hep-ph/9705372]; 

[24] I. General, D. Gomez Dumm and N. N. Scoccola, Chiral phase transition in a covariant 
nonlocal NJL model, Phys. Lett. B 506, 267 (2001) [hep-ph/0010034]; 

[25] T. Sakai and S. Sugimoto, Low energy hadron physics in holographic QCD, Prog. Theor. 
Phys. 113, 843 (2005) [hcp-th/0412141]. 

[26] T. Sakai and S. Sugimoto, More on a holographic dual of QCD, Prog. Theor. Phys. 114, 1083 
(2005) [hcp-th/0507073]. 

[27] E. Antonyan, J. A. Harvey, S. Jensen and D. Kutasov, NJL and QCD from string theory, 
arXiv:hcp-th/0604017. 

[28] J. A. Harvey, Chiral symmetry breaking and intersecting D-brane systems, Nucl. Phys. Proc. 
Suppl. 171, 243 (2007). 

[29] A. Dhar and P. Nag, Intersecting branes and Nambu— Jona-Lasinio model, Phys. Rev. D 79, 
125013 (2009) [arXiv:0901.4942 [hep-th]]. 

[30] J. F. Donoghue, E. Golowich and B. R. Holstein, Dynamics of the standard model, 
Cambridge Universe Press (1992). 

[31] R. Jackiw, Functional evaluation of the effective potential, Phys. Rev. D 9, 1686 (1974). 

[32] H. Kleinert, Particles and Quantum Fields (chapter 15), 

http:/ /users. physik.fu-berlin.de/ kleinert/b6/psfiles/Chapter-15-direffac.pdf. 



- 23 - 



